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ABSTRACT 

We combine H i 21cin observations of the Milky Way, M31 , and the local galaxy population with QSO 
absorption-line measurements to geometrically model the three-dimensional distribution of infalling 
neutral gas clouds ( "high- velocity clouds", HVCs) in the extended halos of low-redshift galaxies. We 
demonstrate that the observed distribution of HVCs around the Milky Way and M31 can be modeled 
by a radial exponential decline of the mean H i volume filling factor in their halos. Our model suggests 
a characteristic radial extent of HVCs of i?haio ~ 50 kpc, a total Hi mass in HVCs of ^ 10^ Af©, and 
a neutral-gas-accretion rate of ^ 0.7Moyr~^ for M31/Milky- Way- type galaxies. Using a Holmberg- 
like luminosity scaling of the halo size of galaxies we estimate i?haio ~ 110 kpc for the most massive 
galaxies. The total absorption-cross section of HVCs at z « most likely is dominated by galaxies with 
total Hi masses between 10*-^ and 10^° Mq. Our model indicates that the Hi disks of galaxies and 
their surrounding HVC population can account for 30 — 100 percent of intervening QSO absorption-line 
systems with log A^(Hi) > 17.5 at z « 0. We estimate that the neutral-gas accretion rate density of 
galaxies at low redshift from infalling HVCs is dMm/dt/dV ~ 0.022 Afoyr"^ Mpc~'^, which is close 
to the measured star-formation rate density in the local Universe. HVCs thus may play an important 
role for the on-going formation and evolution of galaxies. 
Subject headings: ISM: clouds - quasars: absorption lines - Galaxy: halo 



1. INTRODUCTION 

One crucial aspect of galaxy formation and evolution 
concerns the continuous infall of intergalactic gas onto 
galaxies. While it is clear that galaxies do accrete sub- 
stantial amounts of gas from intergalactic space to power 
star formation, the exact way of how galaxies get their 
gas is still a matter of debate. In the conventional sketch 
of galaxy formation and evolution gas is falling into a 
dark matter (DM) halo and then is shock-heated to ap- 
proximately the halo virial temperature (a few 10^ K, 
typically) , residing in quasi- hydrostatic equilibrium with 
the DM potential well (Rees & Ostriker 1977). The gas 
then cools slowly through radiation, condenses and set- 
tles into the center of the potential where it forms stars as 
part of a galaxy ('hot mode' of gas accretion). It has been 
argued, however, that for smaller DM potential wells the 
infalling gas may radiate its acquired potential energy 
at much lower temperatures (< 10^'^ K, typically), so 
that one speaks of the 'cold mode' of gas accretion (e.g.. 
White & Rees 1978). For the cold mode of gas accre- 
tion the star-formation rate of the central galaxy is di- 
rectly coupled to its gas-accretion rate (White & Frenk 
1991). Numerical simulations indicate that for individ- 
ual galaxies the dominating gas-accretion mode depends 
on the mass and the redshift (e.g., Birnboim & Dekel 
2003; Keres et al.2005). The general trend for z « is 
that the hot mode of gas accretion dominates for mas- 
sive galaxies with DM- halo masses > 10^^ Mq, while the 
cold accretion mode dominates for galaxies with smaller 
DM-halo masses (e.g., vandeVoort et al. 2011). 

Independently of the theoretically expected gas- 
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accretion mode of galaxies it is known since a long time 
that galaxies at low and high z are surrounded by large 
amounts of neutral and ionized gas that partly originates 
in the IGM. This material is complemented by neutral 
and ionized gas that is expelled from the galaxies as part 
of galactic fountains, galactic winds, and from merger 
processes (see, e.g., Richter 2006 for a review). Because 
the interplay between these circumgalactic gas compo- 
nents is manifold and the gas physics of such a turbu- 
lent multi-phase medium is complex, the circulation of 
neutral and ionized gas in the inner and outer halos cur- 
rently cannot be modeled in full detail in hydrodynam- 
ical simulations. To improve current models of galaxy- 
evolution models it is of imminent importance to quan- 
tify the amount of cool, neutral gas in and around galax- 
ies from observations and search for observational strate- 
gies to separate metal-deficient infalling intergalactic gas 
from metal-enriched gaseous material that is circulating 
in the circumgalactic environment of galaxies as a result 
of fountain processes and galaxy mergers. 

In the Milky Way and other, very nearby spiral galax- 
ies (e.g., M31), the infall of neutral gas onto to the disks 
can be observed directly by Hi 21cm observations of 
extraplanar gas clouds that move through the halos of 
these galaxies. For the Milky Way, the so-called "high- 
velocity clouds" (HVCs) represent the prime candidates 
for neutral gas that is being accreted onto the Milky 
Way disk. HVCs represent high-latitude gaseous struc- 
tures (located in the Galactic halo) observed in Hi 21cm 
emission at high radial velocities, |wlsr| > 100 kms^^ 
(e.g.; Wakker & van Woerden 1998; Richter 2006). Halo 
clouds with somewhat smaller radial velocities in the 
range |wlsr| = 50 — 100 kms~^ are commonly referred 
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to as "intermediate- velocity clouds" (IVCs). Through- 
out this paper, we will use the expression "HVC" for 
all neutral halo clouds (including IVCs) , if not otherwise 
stated. The total neutral gas mass of the Milky Way's 
HVC population is on the order of 10^ Mq and the to- 
tal accretion rate of neutral gas in the form of HVCs 
has been estimated to be ^ 0.5 Mq yr^^ (eg-, Wakker et 
al. 2007, 2008; Wakker 2004). Also M31 exhibits a popu- 
lation of neutral halo clouds in a similar mass range (for 
simplicity, hereafter also referred to as "HVCs" ; Thilker 
et al. 2004). These observations, together with Hi 21cm 
measurements of other nearby galaxies (e.g., Sancisi et 
al. 2007), imply that HVCs represent a common phe- 
nomenon in the local Universe. 

Another important method to study the gaseous out- 
skirts of galaxies and their relation to the cosmic web is 
the analysis of intervening absorption-line systems in op- 
tical und uhraviolet (UV) spectra of QSOs (AGN). QSO 
absorption spectroscopy allows us to detect both neutral 
and ionized gas in the intergalactic medium and the ha- 
los of galaxies over eight orders of magnitude in column 
density and over more than 90 percent of the age of the 
Universe. It is therefore a particularly sensitive method 
to explore the multi-phase nature of circumgalactic gas 
and its origin (e.g., Bergeron & Boisse 1992; Steidel 1995; 
Charlton & Churchill 1998; Ding et al. 2005; Richter et 
al. 2011). 

In this paper, we combine H i 21cm observations of the 
HVC population of Milky Way and M31 with 21cm data 
of the local galaxy population and QSO absorption-line 
measurements to study the three-dimensional distribu- 
tion of (partly) neutral gas structures in the extended 
halos of low-redshift galaxies. While several previous 
studies have linked the Galactic HVC population to in- 
tervening QSO absorbers based on various arguments 
(e.g.. Blitz et al. 1999; Charlton, Churchill & Rigby 2000; 
Pisano et al. 2004; Mshar et al. 2007; Schaye, Carsweh 
& Kim 2007; Narayanan et al. 2008; Richter et al. 2009; 
Stocke, Keeney & Danforth 2010; Ribaudo et al. 2011), 
a detailed geometrical HVC model that connects key ob- 
servables (column densities, covering fractions, absorber 
number densities) of HVCs and intervening absorbers 
with the gas-accretion rates of galaxies at low redshift 
has not been presented so far. The main goal of the 
present study is to model the radial distribution of HVC 
analogs in halos as a function of galaxy mass and size, 
determine their absorption-cross section at z w 0, and 
estimate the neutral-gas accretion rate of galaxies at low 
redshift. 

We refrain from including Hi 21cm observations of 
HVCs from other galaxies beyond the Local Group in 
our study, because these observations are strongly lim- 
ited in sensitivity for detecting diffuse neutral halo gas 
(e.g., Oosterloo et al. 2007). In addition, beam-smearing 
effects are known to smoothen out the true spatial dis- 
tribution of individual halo clouds in more distant galax- 
ies, so that the observed filling factor of Hi 21cm in 
these systems does not provide a realistic estimate for 
the absorption-cross section of this gas. 

The paper is organized as follows: in §2 we discuss Hi 
21cm observations of the HVC population of the Milky 
Way and M31. In §3 we develop a simple model for the 
radial distribution of neutral gas in the halos of these 
galaxies. In §4 we generalize our HVC model for the local 



TABLE 1 

Properties" of prominent Milky 

\\AY II\T'S 



HVC name 


J HVC 






[kpc] 


Complex A 


0.01 


8-10 


Complex C 


0.04 


~ 10 


Complex H 


0.01 


> 5 


Magellanic Stream 


0.04 


~ 50 


Complexes WA, WB 


0.01 


8-20 



References: Wakker (2001,2004); 
Wakker et al. (1999, 2007); B.P. Wakker 
(priv. comm.); Thom et al. (2006, 2008); 
Gardiner & Noguchi (1996) 
^ Sky-covering fraction of HVC gas 
Distance 

galaxy population, based on information on the H i mass 
function of low-redshift galaxies and the absorption-cross 
section of Damped Lyman a absorbers (DLAs) and Ly- 
man limit systems (LLS) at ^; ~ 0. The relation between 
HVCs and intervening QSO absorbers is discussed in §5. 
In §6 we present the conclusions from our study. 

2. HIGH- VELOCITY CLOUDS IN THE HALOS OF THE 
MILKY WAY AND M31 

2.1. HVCs in the Milky Way 

The HVC population of the Milky Way has a total 
sky covering fraction of /c « 0.30 for column densities 
iV(Hi) > 7 X 10^'^ cm-2 and /c w 0.15 for column densi- 
ties A/' (Hi) > 2 X 10^^ cm-2 (Wakker 2004 and references 
therein). The largest Milky HVC is Complex C, which 
covers ^ 1500 square-degree on the sky (/c ~ 0.04). 
With a low metallicity of ^ 0.15 solar (e.g., Fox et 
al. 2004; Richter et al. 2001) and a distance of £» ~ 10 
kpc (Wakker et al.2007; Thom et al.2008). Complex C 
most likely is a cloud that is being accreted from the 
IGM or from a satellite galaxy. The Magellanic Stream 
(MS) also covers an area of ~ 1500 square-degree, but 
most likely has a distance as large as D ~ 50 kpc (Gar- 
diner & Noguchi 1996). The MS represents a tidal fea- 
ture expelled from the Magellanic Clouds as they move 
through the extended Milky Way halo. With its large dis- 
tance and its stream-like shape the MS is clearly distinct 
from most of the other Galactic HVCs, which predomi- 
nantly are less extended and located at distances < 15 
kpc from the disk (Wakker et al. 1999, 2007, 2008; Thom 
et al.2006, 2008). 

Other prominent Galactic HVCs are Complex A, 
Complex H, the Anti-Center Cloud, and Complexes 
WA— WE. Covering fractions and distances for some of 
these complexes (as far as known) arc summarized in 
Table 1. Because accurate distance information on the 
Milky Way HVCs is still limited, the total Hi mass in 
HVCs and the Hi mass-accretion rate ist not well con- 
strained. The currently available information implies 
Mm HVC 3 X 10^ Mq and dMni Yi\c/dt - 0.5 Mq yr"! 
(Wakker et al.2007; 2008; Wakker 2004). The Hi col- 
umn densities in HVCs follow a column-density distribu- 
tion function (CDDF) in the form f{Nm) oc N^f with 
j3 = 1.42 for log N{lli) > 18 (Lockman et al. 2002). 
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Fig. 1. — Left: distribution of HVCs (green) around M31, adopted from 21cm data of Thilker et al. (2004). The optical disk of M31 
is indicated with the black ellipse, the H I disk is shown as gray-shaded ellipse. Note that the shown HVC distribution is incomplete for 
small radii (r < 15 kpc) due to confusion of Hi halo gas with the Hi disk of M31. Right: projected covering fraction of HVC gas around 
the M31 disk, /hvCi a- function of the projected radius, r, based on the HVC distribution shown in the left panel. The covering fraction 
can be fitted by an exponential in the form 2.1 exp(— r/12), as indicated with the red solid line. 

2.2. The M31 HVC population 

As has been shown by Braun ct al. (2009), M31 has an 
H I disk that extends to radii of r « 30 kpc for column 
densities above log A^(Hi) > 20.3. Beyond r « 30 kpc 
the Hi disk appears to be truncated (Braun et al. 2009; 
Braun & Thilker 2004). The HVC population of M31 was 
studied in detail with the Green Bank Telescope (GBT) 
by Thilker et al. (2004) . These authors detected more 
than 20 individual HVCs around M31 with Hi column 
densities > 5 x 10^^ cm~^ and estimated a total Hi mass 
of the HVCs of 3 - 4 X 10'' Mq (see also Thilker, Braun 
& Westmeier 2005). 

In the left panel of Fig. 1 we show the distribution of 
HVCs around the M31 disk, based on the GBT Hi 21cm 
contour map presented by Thilker et al. (2004). The 
HVC population of M31 reaches out until ^ 50 kpc with a 
strongly decreasing HVC covering fraction towards larger 
radii. We assume D = 785 ± 25 kpc as distance for M31 
(McConnachie et al. 2005), so that absolute distance esti- 
mates derived from angular coordinates are uncertain by 
~ 3 percent. Some interesting conclusions about the ra- 
dial distribution of neutral gas around M31 can be drawn 
from this HVC distribution map. Fig. 1, right panel, 
shows the radius-dependent projected covering fraction, 
/hvc(?'), of the HVC population of M31 plotted against 
r (indicated by the filled boxes), where r is the projected 
radius. To calculate /hvc(^) we have resampled the M31 
HVC map of Thilker et al. (2004) and have transformed 
the (a, 5) coordinate system into a polar coordinate sys- 
tem with coordinates r and centered on M31. For each 
ring with radius r and thickness r ± Ar the parameter 
/hvc(''') then was derived by comparing the area cov- 
ered by HVC gas with the total ring area. The error 
bars for /hvc('') shown in Fig. 1 have been calculated as- 
suming Poisson-like statistics. Starting from /hvc ~ 0.5 
at r = 15 kpc the covering fraction decreases to values 
less than 0.05 for radii larger than r — 45 kpc. This 
trend for /hvc ('') can be fitted by an exponential in the 
form /hvc('') — xexp{—r/y) with x — 2.1 ± 0.2 and 
y = 12.0lQ'g, as shown by the solid red line in the right 



panel of Fig. 1. An exponential fit to the extraplanar Hi 
features of M31 was also favored by Braun & Thilker 
(2004), who analyzed lower-resolution 21cm data of M31 
from the Westerbork Synthesis Radio Telescope (WSRT). 

For the inner regions of the M31 halo at r < 15 kpc 
no HVC data are available. As discussed by Thilker et 
al. (2004), this does not imply that this region is devoid 
of HVC material. The lack of data for the disk-halo in- 
terface region at r < 15 kpc (Fig. 1, left panel) rather 
indicates the incompleteness of the HVC map for small 
radii because of the confusion of neutral halo gas with 
the Hi disk of M31 together with the stringent selec- 
tion criteria defined by Thilker et al. to unambiguously 
identify HVC features. In the Milky Way, the disk-halo 
interface at r < 15 kpc is filled with large amounts of neu- 
tral gas that gives rise to 21cm emission at intermediate 
and high velocities (e.g., Wakker 2004). Nearby edge- 
on galaxies such as NGC891 also exhibit large amounts 
of neutral gas in the disk-halo interface region extend- 
ing several kpc above and below the disk (see Sancisi et 
al. 2008). 21cm measurements of NGC 891 show that the 
(projected) covering fraction of neutral gas is /hi = 1 
for vertical distances d < 10 kpc to the midplane of the 
NGC 891 disk (Oosterloo et al. 2007). If we extrapolate 
/hvc for M31 to small radii using the exponential defined 
above, /hvc becomes unity for r < 9 kpc, in line with 
the extraplanar gas distribution observed in NGC 891. If 
we define as the radius beyond which the projected 
covering fraction of HVCs falls below 3 percent, we de- 
rive for M31 a value of ra = 50 ± 6 kpc. The la error 
refiects the uncertainties in the exponential parameters 
for favcir)- 

As shown by Braun & Thilker (2004), the high- 
resolution 21cm HVC data of M31 from Thilker et 
al. (2004) follows a "standard" H i CDDF in the form 

f{Nm) oc Nj^f with /3 w 1.5 in the column density range 
log 7V(H i) = 18 — 20. The slope is very similar to the one 
derived for the Milky Way HVCs (Lockman et al. 2002; 
previous subsection) and is in good agreement with val- 
ues derived for low-redshift QSO H i absorption-line sys- 
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TABLE 2 

Model parameters for Milky Way/M31 HVCs 



Parameter 


Description 


Value 


Input 


ftHVC 


scale height of HVC population 


6.67 kpc 


/v,0 


central HVC volume fraction 


0.0185 


(nm) 


mean H i volume density 


0.1 cm-3 


^IVC 


upper z-height limit for IVCs 


5 kpc 


^'infall 


infall velocity 


100 kms-i for HVCs 






50 kms-i for IVCs 



Output 

/hvC ('') projected HVC covering fraction for radii r > 9 kpc 2.1 exp(— r/12) 

/hvc(^) projected HVC covering fraction for radii r < 9 kpc 1 

(/ivcMw) mean IVC covering fraction from inside the sphere 0.30 

(/llvc;.Mw) mean HVC covering fraction from inside the sphere 0.29 

(/hvc) mean HVC/IVC covering fraction from an outside vantage point 0.21 

r3 radius for which (/hvc) > 0.03 50 kpc 

(A'^(Hl)) mean Hi column density in HVCs/IVCs from inside the sphere 1.3 X 10^'^ cm~' 

.^Hi.tot total neutral gas mass in HVCs/IVCs at r < r-^ 1.2 x 10*^ Mq 

dM^ii^tot / dt total neutral gas mass accretion rate for gas at r < ra 0.74M© yr"'^ 



terns (see Sect. 5.1). 

3. GEOMETRICAL MODELING OF HVCS IN THE LOCAL 

GROUP 

3.1. Modeling setup 

To combine the observational information on the HVC 
population of the Milky Way and M31 we have developed 
the custom-written numerical code halopath, which is 
based on a simple geometrical model assuming spheri- 
cal symmetry. The code allows us to model the radial 
distribution of gas in the halos of galaxies, its mass dis- 
tribution and its absorption-cross section from any given 
vantage point inside and outside the sphere. Because of 
the unknown size distribution of HVCs we here do no 
attempt to model individual H i clouds as HVC analogs, 
but instead consider the volume-filling factor of neutral 
gas as main input parameter, from which all relevant 
physical quantities (e.g., projected covering fraction, Hi 
column density, Hi mass) can be easily obtained and 
compared to observations. 

To support the model with the necessary observational 
data, we assume that the HVC populations of the Milky 
Way and M31 are identical in a statistical sense (same 
radial distribution and same volume filling factor of gas 
with log 7V(Hi) > 17.5). In view of the similarity of 
both galaxies in terms of morphology, mass, luminosity, 
etc. this assumption is justified. As mentioned above, 
the key parameter that describes the spatial distribution 
of HVCs in our spherical model is the radius-dependent 
volume-filling factor of optically thick HVC gas with log 
iV(Hi) > 17.5, /v,HVc(-R), where R is the physical ra- 
dius (compared to the projected radius r). Because the 
projected (area) covering fraction of HVCs in the M31 
halo can be described by an exponential (see above), we 
assume that /v,HVc(-R) follows an exponential, too, so 
that we write: 

,/v,HVC = ./v,0 CXp (-i?/ft,HVc)- (1) 

In this equation, /v,o is the volume-filling factor in the 
center of the sphere and /ihvc is the scale height of the 
HVC population. Another important parameter in our 
model is the mean Hi volume density (nni) in the HVC 
gas, as (riHi) together with fv,Hvc{R) determines the 



neutral gas mass in HVCs per radial bin {— radial vol- 
ume element), Mui{r) = /v,hvc(-R) ('^hi) MTOh, with toh 
as hydrogen mass and /x as a factor that corrects for 
the presence of helium and heavy elements in the gas. 
In addition, (nni) determines the mean Hi column den- 
sity, (A^(Hi)), measured along any given line of sight 
through the halo, since (A'"(H i)) = (nHi) d, where d is 
the absorption-path length through the halo. The pa- 
rameters (A^(Hi)) and (nni) are constrained by 21cm 
observations and ionization models of Galactic IVCs and 
HVCs. Note that in this study we do not model the gas 
physics in the halo clouds, but consider only the spatial 
distribution of neutral halo gas, its Hi column density 
distribution, its total mass, and its infall rate. 

To model the sky covering fraction of neutral gas from 
a vantage point inside the sphere (i.e., to model the pro- 
jected neutral gas distribution in the Milky Way halo 
from the position of the sun) we introduce additional 
constraints based on results from 21cm observations of 
IVCs and HVCs. First, we only consider H i gas as neu- 
tral halo gas if it is located at vertical distances z > 600 
pc from the midplane of the disk (whose orientation can 
be chosen in the model). Second, we separate neutral 
halo gas close to the disk (small z-hcights) from more 
distant halo clouds via the parameter zivc • This param- 
eter enables us to distinguish between IVCs (z < zivc) 
and HVCs {z > Zivc)- Finally, the parameter i;infaii de- 
fines the infall velocity of the neutral halo gas; Uinfaii 
can be a function of i? or z, or can be chosen to be 
constant for IVCs and HVCs, respectively. The neutral 
gas (mass) accretion rate per radial bin then is given by 
dMni{r)ldt = MHi(r) Vi„fau/r. 

3.2. Modeling results 

Using our geometrical model we are able reproduce 
the observed properties of the HVC population of the 
Milky Way and M31 using appropriate values for the 
above discussed input parameters. Our favorite model is 
summarized in Table 2. In this model, the scale- height 
of the HVC population is ft-nvc = 6.671q;45 kpc and 
the central volume-filling factor of neutral halo gas is 
/v,o = 0.0185 ±0.0036. The errors have been determined 
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numerically; they reflect the Icr error range of the ex- 
ponential parameters for /hvc in M31 (Sect. 2.2). In 
our model we assume that IVCs and HVCs are sepa- 
rated at a ^;-height of zi\c = 5 kpc, in line with what is 
known about the distribution of IVC and HVC distances 
in the Milky Way halo (Wakker et al. 2007, 2008; Thom 
et al.2006, 2008). 

Only these three parameters (/ihvC: /v.Oj and zivc = 
5) are required to reproduce the observed exponential de- 
cline of the projected covering fraction of the HVC pop- 
ulation of M31 (/hvc = 2.1 exp (-r/12) for radii r > 9 
kpc) from a vantage point outside the sphere and the 
observed covering fractions of ^ 30 percent for IVCs 
and HVCs from a vantage point inside the sphere (see 
Sect. 2). Note that for r < 9 kpc the projected covering 
fraction is set to unity. The characteristic radial extent of 
the HVC popiilation around both galaxies is = 50 ± 6 
kpc. The mean projected covering fraction of neutral 
IVC/HVC gas in the halo region is (/hvc) = 0.21 ±0.02. 

OlDservations of H l, metal ions and molecular hydro- 
gen suggest that the neutral hydrogen volume densities 
in IVCs and HVCs may span a large range over at least 
three orders of magnitude (10^^ < nm < 10 cm~'^; e.g., 
Wakker et al. 1999, 2004; Richter et al. 2003a, 2003b, 
2009; Sembach et al. 2001). For our model we adopt 
a value of (nni) = 0.1 cm~^, which is close to the one 
that has been derived for high-velocity cloud Complex C 
(Wakker et al. 1999). We consider this value as a realistic 
estimate for the volume-averaged mean neutral hydrogen 
density in IVCs and HVCs. As infall velocity we adopt 
Winf^u = 100 kms-i for HVCs and 50 kms'^ for IVCs, 
assuming that the (vertical) infall velocities towards the 
disk are comparable to the observed radial velocities of 
IVCs and HVCs. Since the infall velocity of HVCs is de- 
termined by the balance between the gravitational force 
and the ram-pressure force provided by the surrounding 
hot coronal gas (e.g., Benjamin & Danly 1997; Briins & 
Mebold 2004), a more precise modeling of Winfaii would 
also require the modeling of the density distribution of 
hot halo gas, which is beyond the scope of this study. 
However, from grid-based hydrodynamical simulations of 
HVCs Heitsch & Putman (2009) derive infall velocities 
that are very similar to the velocities adopted by us. 

With the above given values for rini and Uinfaii the 
mean Hi column density in IVCs and HVCs (from the 
interior view) is (-/V(Hl)) = 1.3 x 10^^ cm-^, the total 
H I mass in the halo ai r <r^ is MHi,tot = 1.2 x 10* Mq, 
and the total H i mass accretion rate for gas at r < rs is 
dM\iixoi/ dt = 0.74 M0yr^^. These values are in excel- 
lent agreement with the observations (see §2). 

4. MODELING OF GALAXY ABSORBERS 

4.1. Galaxies and their absorption characteristics 

In the previous section we have demonstrated that it 
is possible to reproduce the statistical properties of the 
HVC population of the Milky Way and M31 using a 

model that is based on very simple geometrical assump- 
tions. In the following, we want to generalize our HVC 
model for the local galaxy population to constrain the 
absorption cross section of HVCs in the local Universe. 

Galaxies and their circumgalactic gaseous environment 
can be traced by intervening absorption lines of H l and 
metal ions in the spectra of distant QSOs and AGN. The 
strongest intervening neutral-gas absorbers are the so- 



called Damped Lyman a Absorbers (DLAs), which have 
Hi column densities log iV(Hi) > 20.3. These systems 
contain a substantial fraction of the neutral gas mass in 
the Universe (Wolfe et al. 1995). Although there still is 
no consistent picture about the host galaxies of DLAs, 
observations suggest that a mixed population of galax- 
ies contribute to the absorption-cross section of DLAs at 
2 « (e.g., Turnshek et al. 2001; Chen & Lanzetta 2003; 
Rao et al. 2003). From H i 21cm observations of the local 
galaxy population Zwaan et al. (2005) concludes, how- 
ever, that the total DLA cross section at ^ w is dom- 
inated by the gaseous disks of L* and sub-L* galaxies 
with Hi masses > 10^ Mq. 

Based on Hi 21cm observations of the Milky Way, 
M31, and other nearby galaxies it is expected that 
neutral gas absorbers in the extended halos of galax- 
ies (i.e., HVC analogs) have Hi column densities below 
that of DLAs (e.g., Wakker 2004; Thilker et al. 2004; 
Lockman et al. 2002; Murphy et al. 1995). Using the 
common absorber classification scheme, halo absorbers 
therefore are expected to be seen as so-called sub-DLAs 
(19.0 <log Af(Hi) < 20.3) and Lyman-Limit Systems 
(LLS; 17.2 < log iV(Hi) < 19.0). A large fraction of 
the LLS at low redshift therefore may represent distant 
analogs of the HVCs seen around the Milky Way and 
M31 (Richter et al.2011). 

As for HVCs, the H l column-density distribution func- 
tion of intervening QSO absorbers at low and high red- 
shift below the DLA column-density limit can be fitted 
by a power-law in the form /(A^hi) oc N^f . For QSO 
absorbers, /3 has values between 1 — 2, depending on red- 
shift and the column density interval chosen (see Lehner 
et al. 2007). Unfortunately, /3 is poorly constrained for 
log A''(Hi) > 16 at 2: = due to the limited amount of 
low-redshift Hi absorption-line data in the UV. In con- 
trast, for high z there exists a large data base that allows 
us to constrain /? at a relatively high accuracy (Ribaudo, 
Lehner & Howk 2011). 

The incidence of intervening DLAs, sub-DLAs and 
LLS in QSO spectra, usually expressed by the quantity 
dAfhhs/dz, the number of optically thick Hi absorbers 
per unit redshift, can be obtained from the integration of 
the H I CDDF over the appropriate column-density range 
(in our case log A''(Hl) > 17.5). Moreover, dJ\fi,Ls/dz is 
proportional to the space density of galaxies, rigai, and 
the mean geometrical cross section of optically thick H I, 
{Am), in these galaxies : 



LLS 



dz 



-L 



fiNm)dNm = 



ArHI>JVLLS 



H{z) ■ 



(2) 



We adopt H{z) = {fl^, (1 + z)^ + ^^A)^/^ fl"o = 73 
kms~iMpc"\ = 0.238, and Oa 0.762 (Spcrgel 
et al. 2007). We assume that the covering fraction of 
optically thick Hi gas is (/hvc) = 1 in the (inclined) 
disk of a galaxy and (/hvc) < 1 in the surrounding 
halo, where Hi arises in the form of HVCs. Let Adisk 
be the geometrical cross section of the (inclined) disk 
and Ahaio = 'I'^'^haio ^hc cross section of the (spherical) 
halo region with radius .Rhaio- We then can introduce 
an effective HVC cross section for an individual galaxy, 
^HVCcff = (./HVc)(^haio " ^disk), SO that the total area 
covered by optically thick H i in and around a galaxy is 

^HI = ^disk + ^HVC.eff • 
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Fig. 2. — Left: projected covering fraction of HVCs around galaxies, /hvCi a function of the projected radius, r, for four different Hi 
disk masses (log Mjji disk = 7.75(a), 8.95(b), 9.55(c), 10.15(d); see equation 5). Right: expected number density, dj\f /dz, of optically thick 
Hi absorbers originating in galaxy disks (log Af(Hl) > 20.3; gray-shaded area) and disks plus HVCs (log A'^(Hl) > 17.5; green-shaded area) 
per galaxy Hi mass bin. 

have studied in detail these and other properties of 50 
nearby galaxies in the context of the low-redshift DLA 
population. Rosenberg & Schneider find that the total 
(inclination-corrected) Hi cross section of gas disks in 
their galaxy sample in the DLA column-density range 
(log N > 20.3) is given by Adisk = nab /A, where a 
and b are the major and minor axis parameters (in kpc) 
derived at the DLA column density limit (see also Rosen- 
berg & Schneider 2003, their appendix.) They also find 
a remarkably tight correlation between A^i^]^ and the to- 
tal Hi mass of the disk, Mhi (in solar mass units), as 
log^disk=logMHi-6.82. 

Since the mean value for 6 in a sample of randomly 
inclined gas disks is expected to be 0.637 a, we can write 
for the mean disk area (Adisk) ~ 0.16 7ra^ = 0.5 or 
log (^disk) ~ 2 log a — 0.3. Combining this with the ob- 
served relation between ^disk and A/hi, we obtain a rela- 
tion between the total H i mass and the disk radius at the 
DLA limit in the form log a = 0.5 log A/hi — 3.26. Finally, 
Rosenberg & Schneider (2003) find a (weak) correlation 
between the H i mass and the J-band luminosity in their 
sample, log A/hi = 0.44 log Lj-f 5.27, where Lj is in solar 
luminosity units. 



From Hi 21cm measurements of the local galaxy pop- 
ulation, Zwaan et al. (2005) have derived a DLA number 
density per unit redshift of {dAf/dz)uLA = 0.045 ±0.006, 
a result that is good agreement with previous esti- 
mates (e.g., Rosenberg & Schneider 2003; {dj\f / dz)Yii^p^ = 
0.053 ±0.013). These measurements provide direct infor- 
mation on the individual geometrical cross section (i.e., 
^disk) of the gaseous disks at log 7V(Hi) > 20.3 in low- 
redshift galaxies for a large range of galaxy morpholo- 
gies and luminosities. These H i surveys are, however, 
not sensitive enough and do not provide sufficient spatial 
resolution to identify a possibly existing HVC population 
in the halos of these galaxies. 

Based on the above-given relations, a QSO sightline 
passing through both an Hi disk and HVC gas of a 
galaxy would show a DLA as absorption signature, while 
a sightline passing only through an optically thick HVC 
would exhibit a sub-DLA or LLS. We then can write 
for the total number density of H i absorbers with log 
iV(Hi) > 17.5 that trace gas disks of galaxies and their 
surrounding HVC population: 



dJ\f \ 

) disk-l-HVC 



( dM 
V dz 



disk 



iHVCcff 



DLA 



Adi 



disk 



(3) 



Our goal is to estimate ^Hvc,cff from a generalized ver- 
sion of our HVC model for the local galaxy population 
(for which Adisk is known). This will enable us to esti- 
mate (c?A/'/(iz)disk+HVC and link the HVC population at 
2 = with the H i column-density distribution function 
of low-redshift QSO absorbers. 

4.2. On the absorption cross section of neutral gas disks 

Before we start to investigate the absorption cross- 
section of neutral gas in the halos of galaxies, it is use- 
ful to briefly discuss the relation between the absorption 
cross section, the total Hi mass, and the scale- length of 
neutral gas disks in local galaxies, as derived from H i 
21cm surveys. Using 21cm data from Arecibo and the 
Very Large Array (VLA) Rosenberg & Schneider (2003) 



4.3. On the absorption cross section of HVCs 
surrounding neutral gas disks 

Surveys of the local galaxy population, such as the one 
presented by Rosenberg & Schneider (2003), show that 
galaxies span several orders of magnitude in parameters 
like Hi mass, optical luminosity, and Hi disk size. One 
crucial question that concerns the cross section of neutral 
gas in the halos of these galaxies is, how the size of the 
gaseous halo of a galaxy is related to the above listed 
parameters. 

The standard approach to scale the size of a galaxy's 
gaseous halo with its luminosity is to adopt a Holmberg- 
like luminosity scaling, so that the halo radius is given 
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TABLE 3 

Properties'* of H i absorbing galaxies 



log Mni.disk 


log Adisk 


''3 


(/hvc) logAnvCeff (<iA/'/(i0)disk+Hvc 


logMni.HVC dMui/dt 


[Mhi in Mq] 


[A in kpc^] 


[kpc] 


[A in kpc^] 


[Mhi in Mq] [Mq yr"!] 



6.7-7.0 


0.03 


2 


0.35 


0.41 


7.1 X 10--* 


3.38 


5.9 X 10-5 


7.0 - 7.3 


0.33 


3 


0.23 


0.88 


1.9 X 10-3 


4.21 


3.7 X 10-4 


7.3 - 7.6 


0.63 


4 


0.21 


0.97 


8.4 X 10-4 


4.82 


1.4 X 10-3 


7.6 - 7.9 


0.93 


7 


0.15 


1.51 


5.9 X 10-3 


5.46 


5.2 X 10-3 


7.9-8.2 


1.23 


9 


0.17 


1.73 


6.1 X 10-3 


5.91 


1.2 X 10-2 


8.2 - 8.5 


1.53 


13 


0.18 


2.05 


1.3 X 10-2 


6.39 


3.1 X 10-2 


8.5 - 8.8 


1.83 


19 


0.18 


2.40 


1.7 X 10-2 


6.85 


7.5 X 10-2 


8.8-9.1 


2.13 


27 


0.19 


2.70 


3.4 X 10-2 


7.29 


0.17 


9.1-9.4 


2.43 


39 


0.19 


3.00 


3.8 X 10-2 


7.73 


0.40 


9.4-9.7 


2.73 


55 


0.20 


3.30 


3.9 X 10-2 


8.15 


0.85 


9.7- 10.0 


3.03 


78 


0.20 


3.61 


3.8 X 10-2 


8.56 


1.69 


10.0 - 10.3 


3.33 


110 


0.22 


3.91 


1.8 X 10-2 


8.95 


3.09 



* Explanations: r-s = halo radius beyond which the projected HVC covering fraction is < 3 percent; Adisk =Hl 
disk area; (/hvc) =incan projected HVC covering fraction for halo region; Anvc^.cft = /HVc;(Ahaio ^ Adisk); 
{dj\f /(i2)disk+HVC = number density of disk/halo absorbers with log A''(Hl) > 17.5 (per galaxy Hi mass bin); 
-'^^HljHVC = total neutral gas mass in HVCs per galaxy; dMm/dt = neutral gas mass infall rate per galaxy. 



by 

i?halo(i) = R. (^-^) , (4) 

where S » 0.2 and iJ* « 110 kpc for B-band luminosities, 
as derived from analyses of intervening Mg II absorbers 

(which trace neutral and ionized gas in disks and halos) 
and their relation to galaxies (Steidel 1995; Kacprzak et 
al. 2008). It can be shown that with the above given rela- 
tions between Lj, Mhi, and Adisk a Holmberg-like lumi- 
nosity scaling as given by equation (4) corresponds to a 
linear scaling of the halo radius with the H i disk radius. 
This is because Mhi oc L^-^'^ (Rosenberg & Schneider 

2003) and a oc My(^, which leads to i?haio ^ L^'^'^ oc a, 
i.e., a Holmberg-like luminosity scaling with 5 = 0.22. 
For the following, we therefore assume a Holmberg-like 
luminosity scaling with S = 0.22, but we parametrize the 
halo size over the relation Rha\o{a) = 7^, where 7 > 1 
and a is the H i disk radius for H i column densities above 
the DLA limit. For M31 7 w 1.7 (see Sect. 2.2). 

To characterize the covering fraction of HVCs, /hvc: 
in galaxy halos, we assume that the exponential decline 
of /hvc observed in M31 reflects a general behavior of 
HVCs in galaxies in the local Universe. We then can 
express the projected covering fraction of HVCs around 
galaxies as a function of its H i disk length in the form 

f ( \ i ^ : r < 0.3 a \ 

/HVC irj ~ I 2.1. exp (-2.5 r/a) : r>0.3a/- 

For M31, a = 30 kpc (Braun et al. 2009; see above). In 
Fig. 2, left panel, we show /hvc(?') for four different Hi 
masses. In analogy to what has been discussed for the 
M31 HVC population, we define rs as the halo radius, 
beyond which /hvc falls below the 3-percent level and 
set i?haio = fa ■ The scaling relation between the disk and 
halo radius in our model then comes out to -Rhaio = 1-69 a 
for all halo radii considered (i.e., 7 = 1.69). 

Since log a = 0.5 log Mhi — 3.26 (see above), equation 
(5) allows us to calculate the sizes of neutral gas halos 
of low-redshift galaxies as a function of their H i mass 
and luminosity. It also enables us to predict the cross 
section and number density of sub-DLAs and LLS that 



represent HVC analogs, their radial distribution, and es- 
timate the neutral-gas accretion rate of galaxies in the 
local Universe. 

4.4. Modeling results 

As input for our generalized HVC model we adopt the 
H I mass distribution of low-redshift DLAs from Zwaan 
ct al. (2005), as derived from a high-resolution 21cm sur- 
vey of the local galaxy population. Using 21cm data 
from the WRST Zwaan ct al. (2005) have studied the Hi 
properties of 355 nearby galaxies and their contribution 
to the local DLA population. The expected values of 
(d/V /dz)uLA for the different Hi masses of the galaxies 
in their sample are shown as gray-shaded area in the right 
panel of Fig. 2. Obviously, galaxies with H I masses in the 
range log Mhi = 8.8 — 10.0 dominate the absorption cross 
section of DLAs at z = 0. 

Based on these data, we derive for each galaxy H i mass 
bin the radius and area of the H i disk and the H i halo 
using the above discussed relations. Note that we here 
do not take into account the possibility that the neutral 
gas disks of galaxies significantly extend below the DLA 
column-density limit. Using our generalized HVC model 
we then calculate for each mass bin the radius-dependent 
(projected) HVC covering fraction (/hvc; equation 5), 
the eff'ective HVC cross section (Auvces), and the ex- 
pected number density of disk+halo H i absorbers (equa- 
tion 3). Moreover, our model calculates for each galaxy 
H I mass bin the mean projected HVC covering fraction 
((/hvc)) J tlie total neutral gas mass in HVCs (Mhi,hvc), 
and the neutral gas mass accretion rate {dMui/dt). For 
the infall velocities we adopt the values for Uinfaii for IVCs 
and HVCs listed in Table 2. The expected number den- 
sities ((iA/'/(iz)disk+HVC as a function of the galaxy Hi 
mass are indicated with the green-shaded area in the 
right panel of Fig. 2. All results are summarized in Table 
3. Based on these results we derive the following relation 
between the Hi disk mass of galaxies, log Mhi (in solar 
units), and the radius of the neutral gas halo, i?haio (in 
[kpc]): 

log iJhaio = 0.5 log Mhi - 3.03. (6) 
If we integrate the values of {dJ\f /rf^;)disk-i-Hvc listed in 
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Table 3 over the entire mass range, we derive a total num- 
ber density of disk/halo absorbers of (dA/'/(iz)disk+HVC — 
0.212. This value is ~ 5 times larger than the number 
density of DLAs at ^; = (Zwaan et al. 2005), suggest- 
ing that the absorption-cross section of HVCs with log 
-^(H i) > 17.5 exceeds that of DLAs by a factor of ~ 4, on 
average. As for DLAs, the total absorption cross section 
of optically thick H i in disks and halos is dominated by 
galaxies in the mass range logMni = 8.8—10.0. As our 
model indicates, the mean (pro jected) covering fraction 
of HVCs in galaxy halos is small, (/hvc) =0.2, typically. 

4.5. Neutral gas accretion-rate density 

Using the mass accretion rates {dMm/dt) listed in Ta- 
ble 3 together with the Hi mass function of the local 
galaxy population (Zwaan et al. 2005) we can estimate 
the neutral gas accretion-rate density of HVCs (mass ac- 
cretion rate per unit volume) at low redshift. We oh- 
tain dMni/rfi/f^F = 0.022 Mgyr-i Mpc"^. Note that 
this value is calculated under the assumption that all 
of the neutral gas in the halos of galaxies is being ac- 
creted onto their disks, independently of its origin inside 
or outside the host galaxies. The above estimate does not 
include the mass of the ionized gas component of HVCs, 
which may contribute substantially to the total mass of 
multi-phase halo clouds (e.g.. Fox et al. 2010; Winkel et 
al. 2011). Also not included are partly neutral gas frag- 
ments with masses and angular sizes below the detection 
limit of 21cm HVC surveys. Such structures are known 
to exist in the Milky Way halo (Richter et al. 2009), but 
their (total) neutral gas mass most likely is small com- 
pared to the large, extended 21cm HVCs. The role of 
ionized gas is further discussed below. 

The value of 0.022 Mq yr~^ Mpe~^ is remarkably close 
to the star-formation rate density at ^; = (p* = 
0.01 — O.O2M0yr~^ Mpc^^, as derived from ultravio- 
let and infrared observational data (Hopkins & Beacom 
2006). Therefore, cold-gas accretion by HVCs possibly 
plays an important (if not dominating) role in feeding 
galaxies at 2: « with gaseous material to power star 
formation. 

The above given value for the accretion-rate density 
can also be compared with recent estimates of the "cold- 
mode" gas accretion-rate densities at z = from eos- 
mological simulations. Using an SPH code, Keres et 
al. (2009) find dM/dt/dV ~ 0.03 Mq yr"! Mpc'^ for 
cold gas that never exceeded a maximum temperature of 
Tmax = 2.5 X 10"'' K. However, using SPH simulations with 
more realistic gas physics vande Voort et al. (2011) derive 
a much lower "cold-mode" gas accretion-rate density for 
galaxies at z = oidM/dt/dV ~ 0.002 Mq yr'^ Mpc"^, 
while for the superordinate DM halos the cold-mode 
accretion-rate is estimated to be dM/dt/dV ~ 0.01 Mq 
yr~^Mpe~^, thus five-times higher. Therefore, only 20 
percent of the cold gas that enters the DM halo in their 
simulation is actually being accreted as cold gas by the 
central galaxy. Unfortunately, these studies do not pro- 
vide information on the absorption cross section of all 
the cold gas in the DM halos and its radial distribiition 
around the galaxies (independently of whether it is being 
accreted or not), so that a detailed comparison between 
the simulation results and our HVC model is not possible 
at this point. 

It needs to be mentioned that gas that is considered 



as "cold" in the cosmological simulations does not neces- 
sarily end up as Hi high- velocity gas that is detectable 
via 21cm observations. It is expected that a substantial 
fraction of the accreted gas that never was heated up 
to the virial temperature of the host halo remains diffuse 
and "warm", i.e., at low densities (nn < 10~^ cm~^) and 
intermediate temperatures (T = 10^ — 10"'' K). The neu- 
tral gas fraction in such warm gas is expected to be low, 
so that it remains unseen in Hi 21cm emission ("warm" 
mode of gas accretion; Heitsch & Putman 2009; Bland- 
Hawthorn 2008). The existence of such a warm, ionized 
gas component in the halos of galaxies is strongly sup- 
ported by the detection of intermediate- and high-ion 
absorption (e.g., from Sim, Cm, Civ, and Siiv) in the 
halo of the Milky Way (e.g., Fox et al. 2006; Sembach et 
al. 1995, 1999) and in the circumgalactic environments 
of other galaxies (e.g., Ribaudo et al. 2011). Also the 
ionized envelopes of 21em HVC complexes represent sig- 
nificant gas reservoirs that need to be considered for a 
realistic estimate of the total (neutral and ionized) gas 
mass that is being accreted by galaxies (e.g., Winkel et 
al. 2011; Fox et al. 2010). The contribution of the ionized 
gas component to the total gas infall rate is difficult to 
determine, however, as the infall velocity of ionized gas 
may be substantially lower than for neutral gas because 
of hydrodynamical effects that affect the ionized cloud 
envelopes, such as gas stripping, turbulent mixing, and 
heat conduction. The same processes also affect the to- 
tal cloud mass and the neutral gas fraction in HVCs and 
thus influence the H i volume filling factor in these clouds. 
These aspects clearly are best resolved through high- 
resolution hydrodynamical simulations (e.g., Heitsch & 
Putman 2009; Vieser & Hensler 2007). 

5. HVCS AND THEIR RELATION TO INTERVENING QSO 
ABSORBERS 

5.1. HVCs and Hi Lya absorbers 

To estimate the contribution of H i HVCs to the num- 
ber density of optically thick Hi Lya absorbers at low 
z, it is necessary to know the Hi CDDF for log iV(Hl) > 
17.5 at low redshift (equation 2). As mentioned above, 
the CDDF at z = is poorly constrained for this column 
density range owing to the fact that high-column density 
H I absorbers are rare and that the amount of the QSO 
absorption-line data in the ultraviolet is limited. Cor- 
belli & Bandiera (2002) have combined absorption-line 
data for intermediate redshifts from Bandiera & Corbelli 
(2001), low-redshift absorption data from Weymann et 
al. (1998), and 21cm Hi emission-line data from Ryan- 
Weber et al. (2003) to construct the Hi CDDF in the 
range log TV (Hi) « 13 - 21. The CDDF presented in 
Corbelli & Bandiera (2002) can be fitted by power laws 
in the ranges log Af(Hi) < 18 and log iV(Hi) > 20, while 
at log A^(Hi) = 18 — 20 there seems to be a plateau in 
f{Nui)- Since this is exactly the column-density range 
of relevance for HVCs, these data imply that the total 
absorption-cross section of HVCs increases only mildly 
with a decreasing H l cutoff column density. To check 
the validity of this interesting feature in /(A^hi) more 
absorption- line data for the range log A^(Hi) = 18 — 20 
are required. To compare the frequency of sub-DLAs 
and LLS at low redshift with the HVC absorption cross 
section obtained from our model we have combined the 
absorption/emission- line data from Corbelli & Bandiera 
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Fig. 3. — Left: Hi column-density distribution function (CDDF) at 2 0, based on results by Corbelli & Bandiera (2002; solid line), 
Zwaan et al. (2005; solid line) and Leliner et al. (2007; dashed line). Right: cumulative number density of Hi absorbers per unit redsliift 
{{dj\f /c'2)>]v(Hl)) ^ function of the cutoff Hi column density. The gray-shaded area indicates the allowed dj\f /dz range for the CDDF 
based on the data of Corbelli & Bandiera (2002) and Zwaan et al. (2005); the gray-shaded area plus the green-shaded area shows the allowed 
dj\f /dz range for the CDDF extrapolated from the Lehner et al. (2007) data. The red filled circle indicates the value for {dj\f / dz) ^[^jf^^uvc 
derived from our HVC model. 



(2002) and Zwaan et al. (2005) and have constructed a 
"hybrid" Hi CDDF at z w 0, in Fig. 3, left panel, shown 
as solid line. 

Independent constraints for /(A^hi) at low redshift 
come from the study by Lehner et al. (2007), how have 
analyzed low- redshift Lya absorbers based on high- 
resolution RST /STIS data. While their study is lim- 
ited to the column density range log 7V(Hi) « 13 — 16.5, 
it is possible to extrapolate their CDDF from the three 
data points between log A^(Hi) = 14.5 — 16.5 (Lehner et 
al. 2007; their Fig. 14) to log (H l) = 20 with a power 
law f{N) oc N^^, where /3 w 1.3 (Fig. 3, left panel, 
dashed line). This extrapolation nicely connects to the 
21cm data from Zwaan et al. (2005), but lies substan- 
tially above the Hi CDDF based on the Corbelli et al. 
results. For the following, we consider both representa- 
tions of the H I CDDF shown in Fig. 3 as plausible input 
parameters to estimate the absorption cross section of 
Hi in sub-DLAs and LLS at 2; w 0. In a recent study, 
Ribaudo, Lehner & Howk (2011) have compiled a large 
sample of LLS, most of them located at high redshifts 
(z < 2.6). Interestingly, their data points for f{N) fit 
well to the Hi CDDF extrapolated from the Lehner et 
al. (2007) data. 

Fig. 3, right panel, shows the (cumulative) number den- 
sity of Hi absorbers per unit redshift, [dAf / dz)y]^(^Ya)-, 
as a function of the cutoff Hi column density, iV(Hi), 
dervied from integrating the two representations of the 
Hi CDDF over all column densities larger than iV(Hi) 
(equation 2). This plot now can be directly compared 
to the results from our HVC model. The gray-shaded 
area indicates the allowed dAf /dz range (0 — 100 percent 
contribution of HVCs to the H i absorber population) for 
the hybrid Hi CDDF (Corbelh & Bandiera 2002; Zwaan 
et al. 2005; left panel, solid line), while the gray-shaded 
plus green-shaded area shows the allowed dM /dz range 
based on the H i CDDF extrapolated from the Lehner et 
al. (2007) data (left panel, dashed hue). The red filled 



circle indicates the value for ((iA/'/(iz)disk+HVC derived 
from our HVC model. Thus, our HVC model predicts a 
value for dN /dz from Hi disk and halo absorbers that 
lies slightly above the value expected from the hybrid 
CDDF, but lies at a ~ 30 percent level of dN /dz pre- 
dicted by the CDDF extrapolated from the Lehner et 
al. (2007) data. We conclude that HVCs in our model 
contribute with 30 — 100 percent to the population of Hi 
Lya absorbers with log iV(Hi) > 17.5 at z w 0. Most 
likely, the contribution of HVCs is clearly less than 100 
percent, as extended neutral and partly-ionized gas disks 
with log A^(H i) < 20.3 and galaxy outflows (as traced by 
strong Mgii absorption; see next subsection) contribute 
to the population of Ly a absorbers in the LLS and sub- 
DLA column density range. 

A more detailed comparison between the absorption 
cross section of HVCs and the local Hi Lya absorber 
population has to await a more precise determination of 
the Hi CDDF at z « from BST / COS data. 

5.2. HVCs and intervening Mgii absorbers 

The Mgll resonant doublet near 2800 A is commonly 
used to study neutral and ionized gas in the outskirts of 
galaxies at 0.3 < z < 2.2 (e.g., Bergeron & Boisse 1991; 
Charhon & Churchih 1998; Ding et al.2005). The so- 
called "strong" Mg II systems represent intervening metal 
absorbers in QSO spectra that have rest-frame equivalent 
widths > 0.3 A in the Mgii A2796 line. They are 
commonly found within 35/i~^ kpc of luminous galaxies 
and thus most likely are related to neutral and ionized 
gas in the disks and halos of low-redshift galaxies. The 
number density of strong intervening Mgii systems in 
the local Universe is expected to be (d7\/'/dz)Mgii ~ 0.5, 
as estimated from extrapolating the redshift evolution 
of strong Mgii systems in the SDSS data from z > 0.3 
down to z w (Nestor, Turnshek & Rao 2005). 

The value for [dN /dz)Mgii \s ^ 2 — 3 times higher 
than (dA/'/dz)disk+HVC estimated from our HVC model. 
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suggesting a substantially larger absorption cross section 
of strong Mg II systems compared to HVC H i absorbers 
with log A^(Hi) > 17.5. This is expected, because with 
an ionization potential of 15 eV Mgii traces both neu- 
tral and ionized gas in the disks and halos of galaxies. 
Kacprzak et al. (2008) estimated a Mgii covering frac- 
tion of (/Mgii) ~ 0.5 for galaxies and their gaseous halos 
from a sample of 37 Mgii selected galaxies at intermedi- 
ate redshift. This is roughly two times the mean value 
for (/disk+HVc) in our HVC model. 

In the Milky Way halo, the Mgll absorption cross 
section in HVCs has not been determined yet, mostly 
because of the lack of appropriate high-resolution NUV 
data. The covering fractions of other low and interme- 
diate ions with strong transitions in the FUV and with 
ionization potentials comparable to that of Mgii (e.g., 
C II, Si ii) arc larger than the covering fraction of opti- 
cally thick H I (Richter et al. 2009), but not large enough 
to explain the high value for {dAf /dz)MgU in the local 
Universe. Recent studies suggest, indeed, that many 
strong Mg II absorbers at low redshift arise in bipolar out- 
flows and galactic winds (e.g.. Bond et al. 2001; Bouchc 
et al. 2011) and thus are related to halo environments 
with larger Mgii covering fractions in galaxies that are 
more actively star-forming than the Milky Way and M31. 
From our estimate for ((iA/'/c?z)disk+Hvc we conclude that 
the contribution of infalling gas clouds (HVCs) to the 
absorption cross section of strong Mg ii absorbers most 
likely is small, but not negligible (< 35 percent). 

5.3. HVCs and intervening Can absorbers 

Next to the Mgii doublet in the near-UV, the two 
Can H&K lines in the optical near 4000 A represent 
valuable tracers for neutral gas in the inner and outer 
regions of galaxies. For instance. Can absorption is fre- 
quently observed in Galactic HVCs (e.g., Richter, West- 
meier & Briins 2005; Ben Bekhti et al. 2008) and is used 
as diagnostic line to derive distance brackets for IVCs 
and HVCs from optical data (e.g., Wakker et al. 2007, 
2008; Thorn et al. 2006, 2008). Based on archival high- 
resolution optical spectra of more than 300 QSOs, ob- 
tained with the VLT/UVES spectrograph Richter et 
al. (2011) have investigated intervening Can absorbers 
at z < 0.5 and their relation to disk and halo gas com- 
ponents. They derive a number density per unit redshift 
of {dJ\f/dz)ca,u = 0.117 ± 0.044 for Can systems with 
log iV(Caii) > 11.65, which is roughly 55 percent of the 
cross section of Hi HVCs with log TV (Hi) > 17.5 esti- 
mated from our model. 

In the Milky Way halo, the covering fraction of Call 
for log 7V(Caii) > 11.65 is ~ 0.2 (Ben Bekhti et al. 2008), 
thus ^ 67 percent of that of 21cm HVCs with log 
iV(Hi) > 17.5. AsdiscussedinRichteretal. (2011), Can 
absorbers with log 7V(Caii) > 11.65 predominantly trace 
neutral gas clouds with log A'^(Hi) > 18.5. This is due to 
the presence of dust in the gas (Ca is strongly depleted 
into dust grains) and the mostly sub-solar metallicities 
of the absorbers. From these numbers it follows that 
intervening Can absorbers arise only in specific regions 
in HVCs, where the Hi column density is large enough 
to compensate for the Ca dust-depletion and abundance 
effects. 

6. CONCLUSIONS 



Several previous studies have discussed a possible link 
between intervening mctal-linc and H I absorbers in QSO 
spectra and the HVC phenomenon in the Milky Way 
(e.g., Blitz et al. 1999; Charlton, Churchill & Rigby 
2000; Mshar et al. 2007; Narayanan et al. 2008; Richter 
et al. 2009; Stocke, Keeney & Danforth 2010; Ribaudo 
et al. 2011). In this study, we move this idea forward by 
bringing together in a quantitative manner the properties 
of HVCs in the Milky Way and M31, QSO absorption- 
line statistics, and neutral gas-accretion rates of galaxies 
in the local Universe. We demonstrate that it is pos- 
sible to explain the star-formation rate density in the 
local Universe through the infall of cold gas in the form 
of HVCs by using a simple geometrical model. In this 
model we project the observed statistical properties of 
the Hi HVCs in the Local Group onto the local galaxy 
population, assuming a Holmberg-like luminosity scaling 
of the halo size. The main results from our modeling are 
summarized in Tables 2 and 3. 

We emphasize at this point that we have not tuned any 
of the input parameters to fit the HVC cross section to 
the H I absorber population or any other above discussed 
QSO absorber observables. Our generalized HVC model 
is based solely on observed properties of the HVC popula- 
tion in the Local Group and the local galaxy population 
without any further assumptions. Note that the exact 
radial dependence of projected HVC covering fraction, 
/hvc('') (here assumed to be an exponential, based on 
the M31 data), is unimportant for the total absorption 
cross section of HVCs, as only the effective area ^HVC,eff 
(i.e., the product of total halo area and the mean HVC 
covering fraction) is relevant for (dAf / dz)iiYc in a sample 
of randomly distributed QSO sightlines. However, be- 
cause the shape of /hvc (*") reflects the volume-filling fac- 
tor and thus the mass distribution of neutral gas around 
galaxies, /hvc('') is closely related to the Hi accretion 
rate of galaxies dMui{r)/dt = Mmir) Vinfaii/r. Thus, 
for a more precise estimate of the neutral-gas accretion 
rate of galaxies at z = it will be important to constrain 
the radial distribution of HVCs around other low-redshift 
galaxies beyond the Local Group using sensitive, high- 
resolution 21cm observations. Such observations need to 
be accompanied by a more precise estimate of the dy- 
namics of neutral gas structures along their infall path 
through the hot coronal gas in the halos of galaxies of 
different masses, e.g., from hydrodynamical simulations 
(Kaufmann et al.2009; Keres et al.2010). 

Since most of the gaseous material that is being ac- 
creted by galaxies at z = may be diffuse, ionized gas 
("warm accretion"; Bland- Hawthorn 2008) rather than 
cold, neutral gas in the form of 21cm HVCs, it will be 
important for future studies to investigate the distribu- 
tion and mass of ionized gas in galaxy halos. Diffuse ion- 
ized gas structures in the extended halos of galaxies are 
expected to have temperatures T < 3 x 10^ K, relatively 
low gas densities, and low neutral gas fractions (i.e., they 
remain unseen in 21cm HVC surveys). Such structures 
can be observed best in absorption in the FUV in the 
lines of low and intermediate ions such as C ii, C ill. Si ii, 
and Sim (Fox et al.2006; Ribaudo et al.2011). Over 
the next few years large amounts of such spectral data 
hopefully will become available from the many ongoing 
observational campaigns with HST/ COS. Results from 
these observations can be easily implemented in our halo 
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model to predict the ionized-gas accretion rate for low- 
redshift galaxies. 

In conclusion, the increasing amount of information 
on the distribution and physical properties of gas in the 
inner and outer halos of galaxies from observations, sim- 
ulations, and semi-analytic models now can be used to 
substantially improve our understanding of gas-accretion 
processes of galaxies in the local Universe. Our study has 
shown one possible way of how to combine information 
from the local HVC population and QSO absorption-line 
systems at z = to investigate these processes for the 
local galaxy population. More detailed studies of this 
kind (e.g., including constraints for ionized halo gas) 
thus could be of great importance to constrain the role 
of gas-accretion processes for the ongoing formation and 
evolution of galaxies at low z and to characterize their 
connection to the cosmic web. 

The author would like to thank Jane Charlton, Glenn 
Kacprzak, and Bart Wakker for helpful comments and 
interesting discussions. 
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